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Abstract 

It may be possible to use operator regularization with Feynman dia- 
grams, which would greatly simplify its use as it has so far been limited 
to the more complicated Schwinger approach. Operator regularization, 
unlike (^-function regularization, is not limited to one-loop order, and pre- 
serves supersymmetry unlike dimensional regularization. In practice the 
use of operator regularization in the context of Feynman diagrams is found 
not to complicate the calculation. 

1 Introduction 

Analytic regularization of quantum field-theory ^ , [2] is not new, but the oper- 
ator regularization approach [3| has not in general been used in conjunction 
with Feynman diagrams, although such use has been implied [5]. Opera- 
tor regularization has the distinct advantage that it can be used with formally 
non-renormalizable theories [S], [7] since the divergences are not reabsorbed, 
but each is removed and replaced by an arbitrary factor; so it might seem well 
worth the effort of having it work with Feynman diagrams, as it has so far been 
limited to the more complicated Schwinger approach. 

While operator regularization does not cure the non-predictability problem 
of non-renormalizability, it has the advantage that the initial Lagrangian need 
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not be extended, as would be the case for regularization approaches that do not 
remove the divergences and so would need the addition of extra terms in the 
original Lagrangian to accommodate them. 

For use with Feynman diagrams, operator regularization in renormalizable 
theories needs to give results equivalent to other methods of regularization; the 
results will not be identical however, since operator regularization, unlike say 
dimensional regularization, removes divergences. The use of operator regular- 
ization with Feynman diagrams to one-loop was tackled previously [5] . 

We begin with a description of operator regularization, how it works, and 
why it should give equivalent results to dimensional regularization. This is 
followed by one and two-loop examples to show how operator regularization in 
practice is no harder then dimensional regularization. 

The use of analytic continuation to deal with the divergences of quantum 
field theory has been criticized [9], but even the formulae use in dimensional 
regularization to deal with all but the logarithmic divergence involve analytic 
continuation. 

2 Operator-Regularization 

The operator regularization scheme is governed by the identity: 

jn / n \ 

= lim — (1 + ai£ + • • • + (1) 

where the a„s are arbitrary, and it is enough that the degree of regularization 
(n) is the loop order. 

There are two separate aspects to this procedure, first the regularization and 
then the continuation, where the divergences are replaced by arbitrary factors; 
these aspects could be separated if so desired. 

2.1 What operator regularization achieves 

Look at operator regularization for a divergent £7^™ 




where the a„s are arbitrary. 

Noting that f2^^^™ has the Laurent expansion 
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n-'-'^ = 4c-„ + • • • + -c_i + Co + 0{e) 

one can now see that the effect of operator regularization is to replace the 
divergent poles by arbitrary constants 

1 

— Q!n 

to yield the finite interpretation 



O " = a„c_„ H h Q!ic_i + Co 



2.2 Generalization 

Operator regularization may be generalized to multiple operators, as appear in 
multi-loop cases. 

A-''... Z-"" = lim — (1 + ai£ + • • • + a„e") — • • • Z"^"™ (2) 

where the a„s are arbitrary, and it is sufficient that n is the loop order; this 
may be more compactly written as 



= A 



-e—k 



2.3 Feynman versus Schwinger 

The operator regularization method was first introduced in the context of the 
Schwinger approach, which while know to be equivalent to the Feynman ap- 
proach, might still leave one asking if operator regularization has the same effect 
in both. So one is lead to asking if operator regularization of the logarithm as 
used in the Schwinger approach 

InQ = - lim — i ^—Q-' ] (3) 

is equivalent to operator regularization as used in the Feynman diagram context, 
namely: 

O-™ = lim — ^-O-^-™ (4) 
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The Schwinger form can be transformed into the Feynman form using 



i-iy 



(m - 1)! dn^^- 



to yield 



n-"' = lim 



d" /(l+e)---(m-l + e)e" 



de' 



(m- 1)! 



which simplifies to 



e^o ds'^ \\ 1/ V m~ljn\ J 



and can be seen to differ from the Feynman form (equation [4]) , so one might 
initially perceive a difference between the results of operator regularization in 
the Schwinger approach from that in the Feynman approach. However, so long 
as one includes all the arbitrary factors in operator regularization 



n-"^ = lim — 



resolutions of zero 



(1 + aieH hane")-rr2" 



V 

they then yield equivalent results. 



(5) 



3 Equivalence at one-loop 

It is one thing that operator regularization should yield the same results for 
the Schwinger and Feynman approaches, but another that operator regulariza- 
tion yield equivalent results as other more well known regularization procedures 
such as dimensional regularization. This has been covered previously in the 
Schwinger approach (10) but the main purpose here is to use operator regular- 
ization with Feynman diagrams. 

Starting from the basic integral that carries the divergences in one-loop 
Feynman diagrams 

r d^'^l 1 _TiA~uj) 1 

J (2^)2- + + 2l.pf ~ r(A) (Af2_p2)^-" 

the other versions follow by differentiating with respect to Py^, so one can con- 
centrate on this one alone. 

One needs to show that operator regularization and dimensional regulariza- 
tion treat the result 
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1 r{A-uj) 1 

(47r)'^ r{A) u^-'^ 
in equivalent ways when divergent, namely in the limit A — lu = 0,-1, 
In operator regularization this is 



(6) 



1 ^.^^ d_ / e{l + as)rie-ny . 



(47r)2 -^^0 de \T{A + e) 
while for dimensional regularization it is 

where n = w — A and having put uj = d/2, where d is the number of space-time 
dimensions. We are interested in comparing these when the original expression 
is divergent, namely when n = 0, 1, . . . etc. 

The operator regularized expression evaluates in this case to (using the help 
of Mathematica |TT] ) 



(-ir 



a-tp{A) + ip{n+l)-\nu] (9) 



(47r)t {A-iy.n\ 

while the dimensionally regularized expression has the Laurent expansion 

^ ' I +\n4TT + ij(n + l) -Inu] +0{e) (10) 



(47r)i (A-iy.n] Ve 
where tp is the polygamma function. 

These agree in form, and so all one-loop results should be equivalent when 
using operator regularization or dimensional regularization. 

3.1 Example: (p^ in 6 dimensions 

This investigation began with the suggcstioij^ to look at the simpler divergent 
integral associated with scalar particles at one-loop in 6 dimensions, and con- 
tinued with probing questions at the two-loop level. 
Looking at this integral (rotated to Euclidean space) 

' ' (11) 



(27r)6 P + m'^(l+ p)2 + 
This integral can be prepared by using the Feynman parameter 'trick 

^private communication with Professor D. G. C. McKeon 
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T(ai+a2) r , - a;)"^-^ 

'ax 



Jo 



to yield 

~ Jo ^ J (27r)6 [^2 + ^2 + p2(l _ ^) + 2Z.p(l - x)f 
One can now proceed by comparing the use of operator regularization for this 
divergent integral to that of dimensional regularization. 

3.1.1 Operator regularization 

Using the generalized operator regularization scheme, the above becomes 



= [ dx lim 4- / 
Jo s^ode J 



d^l / e(l + ae) 



/o e^ode J (27r)6 \[p + w? + p^{l - x) +2l.p{l - x)]'^'^ ^ 
then performing the momentum integrals using the identity 

rf^"/ 1 1 T{A-oj) 



/ 



leads to 



, d eil + ae) Tie - 1) 

dx hm ' 



(13) 



(47r)3 Jo ^"^ ^^-^'o V + t™' + P^^^^ " ^)]'"' / 
where one can now perform the operator regularization limit, using 

d ( e{l + ae) V{e-l) \ . ^. , 

lim — — — = — = u[—a + mu) 

e^o de V r(e + 2) u^-i / 

to yield the finite result 

= J dx (^m^ + p^x{l — x)) (^—a + ln[m^ + p^x{l — x)]) 

This can be evaluated further, but this is not necessary for comparison with 
the result from dimensional regularization, but fj,^ should be included in the 
logarithm, and taken from the arbitrary (a) part, to yield the fixed part 



(4 



1 

/ dx (rri^ + p^x{l — x)) \n 
Jo 



ir + p-j:{l - x) 



(14) 
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and the arbitrary part 



6 



(15) 



(47r)3 

3.1.2 Dimensional regularization 

Regularization might now be done using the dimensional approach; starting 



from the same point (equation 12) with dimensional extension 



= dx 
Jo 



1 



(27r)6-2e + ^2 + p2(l _ ^) + 2/.p(l - x)f 
and again using the identity of equation [13] leads to 

»i 



= r(zi±£) r,, ^ 

Using Mathematica [TT] to help expand about e = 



-u ^ = -T-^ I - - 7 + 1 + ln47r ) M + j—^ + ©(e) 



(47r) 



(47r)3 



where 7 denotes the Euler-Mascheroni constant (0.577. . .); things become, for 
the finite part (again including the jj?) 



leaving the divergent part 



dx (m^ + p^x{l — .t)) In 



+ p'^x{l — x) 



(47r 



^ "'^ - 7 + 1 + In 47r ] ( ' ^ 



(16) 



(17) 



which are seen to agree in form with the result from operator regularization 



(equations 14 and 15) 



3.2 Example: QED 

Applying operator regularization to the three divergent one-loop Feynman dia- 
grams in QED, following Ramond [12j with the Feynman gauge and in Euclidean 
space. 
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3.2.1 One-loop correction to the fermion line 



Starting with the Feynman diagram for the one-loop correction to the fermion 
line (S(p)) (diagrams drawn using JaxoDraw [13] ) 



f7\ 



p-l 

dH 

(2^ ^y-t + m^" I 



Following Ramond [H] this simplifies to 



-ie~ I dx 

'0 



(18) 



(19) 



which is taken as the common starting point for both dimensional and operator 
regularization. 

Now proceeding with operator regularization, following the same general 
route that would be taken with dimensional regularization 



T,{p) = —ie^ / dx lim 



d*l £(1 + a£)7^[j^(l - a;) - m]7^ 



ode J (27r)4 [;2_^^22.+p22.(l_2.)] 



£+2 



using the identity of equation 48 to perform the momentum integrals 



2) [m?x + p'^x{l ~ x)Y 



Mathematica [TT] for the limits 



lim4^ r^(^^l=«-l-ln. 



£^0 de \ r(e + 2) 

and the fact that in 4 dimensions 7^7/^ = —4 and 7;^7p7^ = 27p (from {7^, 7^} 
-2Jp^) yields 



-2i- 



(4^)^ Jo 



dx — x) + 2to) (a — 1 — In [m^x + p^x{l — x)\ ) 



to deliver the fixed part (/i^ taken from the arbitrary a) 
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2i 



(47r)2 



dx 



(1- x) + 2m) In 



rri^x + p'^x{l — x) 



(20) 



and the arbitrary part 



Compare this to the result from dimensional regularization 



(21) 

the finite part 



2i- 



and the divergent part 



dx 



{1- x) + 2m) In 



m?x + p'^x{l — x) 



(47r)2 



--7-1 



I + 4m) —i 



(47r)2 



r2m 



(22) 



(23) 



extra term 

where 7 denotes the Euler-Mascheroni constant (0.577. . .). 

There is a constant difference between these two methods that stems from 
dimensionally continuing the gamma matrices (in dimensional regularization 
alone), but this will be absorbed by the counter terms. 

3.2.2 One-loop correction to the photon line 

Continuing with the diagram for the one-loop correction to the photon line 

i^^Ap)) 

i+p 




Following Ramond [H] this simplifies to 

^^ll^{p) = ^e^iPpiPu ~ 5^uP^) j dx j 



/, + p + m /, + m 



(24) 



dH 



x{l — x) 



(27r)4 [p + m^+p^x{l~x)f 



(25) 
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which is taken as the common starting point for both dimensional and operator 
regularization. 

Proceeding with operator regularization, and again following the same gen- 
eral route that would be taken with dimensional regularization 



d*l e{l + ae) x{l x) 



6^0 de J (27r)4 p + + p^x{l - x)] 



e+2 



Performing the momentum integrals (using the identity of equation 48 1 



—^{Pi^Pi^-Sp^i^p^) [ dx lini f^i^ + ct^ 



r(e) 



again using 



>o de V r(e + 2) w 
yielding the finite part: 



2 {Pfj^Pi' - Sfj.up'^) I dx In 



e^o de V r(e + 2) [to2 + p2^(^i _ j,)]^ 



m + p x{l — x) 



(477) 

and the arbitrary part 



(26) 
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:{Pf^P,^ - S^,t,p'^){a - 1) 



3 (47r)2 

Compare this against the result of dimensional regularization; the finite part 



(27) 



(47r) 

and the divergent part 



2 {Pp-P^^ - Sfj.up'^) I dx In 



m + p x{l — x) 



27r^2 



3 (47r)2 

which are seen to agree in form 



(28) 



(29) 



3.2.3 One-loop correction to the vertex 

Lastly the Feynman diagram for the one- loop correction to the vertex {T p{p, q)) 
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Following Ramond |12j and retaining only the divergent part for this investiga- 
tion, this simplifies to 

where = rn^{x + y) + p'^x{l — x) + q^y{l — y) — 2p.q xy which is taken as 
the common starting point for both dimensional and operator regularization. 

Proceeding with operator regularization, and again following the same gen- 
eral route that would be taken with dimensional regularization 



rW(p,q)--2ze3 



dx 



dy lim 



s^o de J {2tt 



e{l + ae) 



6 + 3 



Performing the momentum integrals (using equation 49 1 



(4^)^ Jo 
then applying 



dx 



^""^ , ,. d fe(l + ae) Ffe) 
dy lim ' 







de V r(£ + 3) (M2 



d fe{l + ae) T{e)\ _ 1 



lim , 

e->o de V r(£ -t- 3) W 



= -(-3 + 2Q;-21nu) 
4 



with ^alplplula = '^lulplp. and 7^7p7p = 27p leads to the finite part 



2ie7p 



dx 



dy In 



M2 



(32) 



where M"^ = m'^{x -\- y) +p'^x{\ — a;) -I- q^y{l — y) — ^p.q xy, and the arbitrary 
part 
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(47r)2 




) 



(33) 



Compare this to the result from dimensional regularization; the finite part 



where again = m^(a; + y) + p^a;(l — a;) + g^y(l — V) — Zp.g xy, and the 
divergent part 



which agree in form, recalling that operator regularization goes further than di- 
mensional regularization in so much as that it actually removes the divergences. 

4 Equivalence at two- loop 

Not surprisingly, demonstrating equivalence at two-loop is somewhat more chal- 
lenging than at one-loop. 

4.1 Example: Two-scoop diagram in 0^ 

Starting with the two-scoop Feynman diagram for 




(34) 




(35) 




In Euclidean-space this diagram 




(36) 



becomes under operator-regularization 




d^q 1 
(27r)4 (g2 + ^2-) 



2 
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It is important to regulate the entire expression and not sub-parts (which would 



give a different result). Evaluating the integrals using the identity of equation 48 
yields 

r(£-l) 1 r(£) 1 



4 (4^)2r(e + 1) (m2)^-i (47r)2r(£ + 2) {m^f 

and finally performing the regularization (with the help of Mathematica [llj ) 
to get the result 

which agrees with the dimensional regularization result from Ramond '12], if 

ai -> - + 1 - 27 
e 

and 



"2^ ^ + -(1-27)^27(1-7) + 



1 1 

e 6 
where 7 denotes the Euler-Mascheroni constant (0.577 . . .). Where we have /i^, 
Ramond has 47r/i2 . So one sees that the divergences have been correctly replaced 
by corresponding arbitrary factors. 

4.2 Example: Setting-Sun diagram in cf)'^ 
Continuing with the first truly two-loop diagram 




Following Ramond [l2] 



S(p) = -y (3m2x(p)+p^if^(p)) (38) 



where 
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K{p) 



dH f d^q 



(27r)4 J (27r)4 (^2 + ^ ^2)[(g _ ; + p)2 + j^2] 



(39) 



and 



KM = 



4.2.1 K{p) 



dH f d^q 



(27r)4 J (27r)4 (^2 + ^2)(;2 _^ ^2) [(g - / + p)2 + ^2 



(40) 



Again it is important to regulate the entire expression and not sub-parts (which 
would give a different result). 



f dH f d^q 



(27r)4 J {2ttY 



[q^ + TV?) {P + m?){{q — I + p)"^ + TV?) 



(41) 



Introducing Feynman parameters 



d'^q 



dH 



(27r)4(g2+^2)2.+2y (2;r)4 

r(2e + 2) 



dx 



r2(£ + 1) 7o ""^ [(/2 + ^2)2. + ((q - ; + p)2 + j„2)(l _ ^)]2e+2 

completing the square in I 



^(2^+2) ^'...^(1-.)^ r 



r2(£ + i) 7o 



1 



d4; 



(2^)4 (g2+^2)2 
1 



2e+2 



(27r)4 [(/ _ + _ ^))2 + ,„2 + (p + g)22.(l _ ^)] 



changing variables T = l — {p + q){l — x) and performing the T integration (using 



the identity of equation 48 1 
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1 r(2£) 

(4^)2 r2(£ + 1)70 



which can be rearranged as 



(27r)4 (q2 + TO2)2e+2 + + ^y^f^l _ j.)] 



2e 



1 r(2e) 



(4^)2r2(e + l) 70 



1 



Now introduce a second Feynman parameter 



2e 



2N2e+2 



(27r)4 

completing the square in q 



[v + q? 



x{l — x) 



y + [q^ + w?){l - y) 



-4e-2 



, \, r(4£ + 2) ^ /-'^ 2e-l(i_ )2.+ l 

(4^)2r2(£ + i)r(2£ + 2) io Jo 



J (2^ 



(2^)4 



(q+py) +py{l-y)+m l-y + 



x(l — x) 



-4£-2 



changing variables g' = q + py and performing the integration (using the 



identity of equation 48 ) 



^At^T ^ f'dx {x{l~x))-' )2e+l 

(4^)4r2(£ + i)r(2£ + 2) 7o Jo 



y 



-I -4e 



x{l — a;) 

The y integral has a divergence at e = which needs to be exposed, using 



15 



and integrating by parts, one finds 



op-rcg 

K{p) 



r(4£) 



(47r)4r2(£ + l)r(2e + 2) 2e Jo 



dx — x)) 



dy 



2e+l 



x(l — x) 



dy y 

-4e 



2e 



(42) 



whicli can be compared to the result from dimensional regularization (Ra- 
mond [H]) 



dim— rcg 

Kip) = 



f dx {x{\-x)r f dyf 



dy 



(4^) 



a;(l — x) 



-2£^ 



(43) 



The finite parts are now much more complicated than in the one-loop cases 
and it is prudent to proceed by showing that the difference between the results 
is zero. Using Mathematica [11] (code given in the appendix) to expand the 
expressions as a Laurent series in e and recalling that operator regularization 
replaces its divergence with a2 and its 1/e divergence with ai one gets, for 
(47r)** times the difference: 



2(1 - 7 + ln47r) Inm'^-I- 



ai 1 

T ^ £ 

ai 012 1 l + 2(ln47r — 7) 
~A ^ ^ 2£2 + 2e 



+ 7^ + ^ + (1 - 27) ln47r + In^ Air 



■0{e) 



One then gets a zero difference, and so agreement, if 



ai -> - + 4(1 - 7 + ln47r) 
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and 



OLi ^ + 



4 8(l-7 + ln47r) 



+8 (1 - 27 + 7^ + 2(1 - 7) ln47r + In^ 47r) + -7r^ 



f £ 

where 7 denotes the Euler-Mascheroni constant (0.577. . .). To agree with Ra- 
mond, we have also left out the /U^ factors. So one again sees that the divergences 
have been correctly replaced by corresponding arbitrary factors. 

4.2.2 K^(p) 

Starting with Kfj_{p) and as before regulating 



KM 



r dH r d^q 
J J2nfj W) 



(p + q- 1), 



{q^ + TV?){P + m^) ((g — I +pY + m^)^ 
Introducing Feynman parameters 



e+l 



(44) 



r d'^q 1 r dH 



{q^ + m'f 

r(3e + 3) /-i^^ a;^(l-x)2^+i(p + g-0^ 



JO 



dx 



r(2e + 2)r(£ + 1) Jo [(/2 + m2)a; + ((9 - / + vf + m2)(l - x)f'^^ 
completing the square in I 



r(2£ + 2)r(£ + i) 

dH 



Y (92 + ^2^+^ 

(p + <? - 



(27r)4 [(z _ (p + 5)(i _ a:))2 + + (p + g)2x(l - x)\ 



3£+3 



letting Z' = Z - (p + g)(l - a;) 



r(3e + 3) 



r(2£ + 2)r(e + l) 



)4 ((Z2+„l2)^+l 

(p + (7)^0; - 



(27r)4 [;,2 _^ ^2 _^ _^ ^)22.(l _ 2;)] 



3e+3 
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performing the I' integral using the identity of equation 48 (dropping odd inte- 
grals) 



1 r(3e + l) 

(47r)2 r(2£ + 2)r(£ + 1) 



rearranging 



1 r(3e + l) 

(47r)2 r(2e + 2)r(£ + 1) 



(27r)4 (^2+ ^2)^+1 

jp + g)p 

[to^ + (p + q)^x{l — x)]^'^~^^ 



dx x~^%l-x)-^ 



d^q 



{P + l), 



{1-x) 



3e+l 



2 I ^2^e+l 



(g2 -|- 77^2^ 



Now introduce a second Feynman parameter 



(27r)4 

completing the square in q 



{(p + 1)^ + 5(f^) y + i<i^ + - y) 



4e+2 



d'^q 

letting q' = q + py 



(P + 



{q + pyY +m?{l-y+ ^(j^) + P^y{l - y) 



4e+2 



(27r)^ 



+m^{l-y+ + J52y(l - y) 



4e+2 
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performing the q' integral using the identity of equation 48 (dropping odd inte- 
grals), one finds 



op— reg 

KM 



r(4£) 



dx X 



(47r)4r(2£ + 2)r2(e + l 

P^2/(l - y) + ( 1 - 2/ + 



-4e 



y 



x{l — x) 



(45) 



which can be compared to the result from dimensional regularization (Ra- 
mond [12]) 



dim— reg 

KM = 



Pp 



r(2£) 

(47r)4-2e 



dx {x{l-x)r' / dyy'{l~y) 



x{l — x) 



-2e 



(46) 



Proceeding as in the K{p) case, by looking at the difference, the two regular- 
izations agree if 



5 + 4(21n2 + ln7r-7) 



where 7 denotes the Euler-Mascheroni constant (0.577. . .). So one again sees 
that the divergences have been correctly replaced by corresponding arbitrary 
factors. 



5 Conclusion 

The above suggests that operator regularization can in fact be used in conjunc- 
tion with Feynman diagrams to all loop orders. The calculation using operator 
regularization is actually somewhat simpler than that using dimensional regu- 
larization, as the gamma matrices are not dimensionally continued when using 
operator regularization. 

While the main purpose of this work is to propose the possibility of using 
operator regularization in the context of Feynman diagrams, it has been noted 
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that the results of operator regularization and dimensional regularization may 
differ [lOj . and that dimensional regularization can have problems respecting 
supersymmetry [14] . 



A Appendix 

A.l Feynman parameters 



D^,'D^' T{ai)ria2)Jo [D,x + 02(1 - x)] 
A. 2 Integrals 

(P'^l 1 1 T{A-iu) 



(48) 

(49) 



(27r)2" + (47r)-r(A) 2 (M2)'4-i-" 

dim— reg op— reg 

A. 3 Mathematica Code for K{p) — K{p) 

This code was composed replacing the operator regularization poles by hand 
(l/e^ and 1/e — > ai) as this made for faster running; it was also found 

that the final double integral ran much faster in an older version of Mathematica 
(version 4), as recent versions are more careful about assumptions. 

arg = 

p-2 y(l-y) + m-2 ( 1-y + y/(x(l-x)) ); 



dimreg = 

-Gamma [2e]/(4Pi)-(4-2e) 1/e * 
(x(l-x))--e y-e D[(l-y) arg-(-2e), y] ; 
opreg = 

-l/(4Pi)-4 Gamma [4e]/ (Gamma [e+1]) "2 l/Gamma[2e+2] l/(2e) * 
(x(l-x))--e y-(2e) D [ (1-y) " (2e+l) arg-(-4e), y] ; 



dimregexp = 

Series [ dimreg, {e,0,0}]; 
opregexp = 
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Series [ opreg, {e,0,0}] ; 

pole2diff = 

l/e"2 Coefficient [ dimregexp, e,-2] - 
a2 Coef f icient [ opregexp, e,-2] ; 
poleldiff = 

1/e Coefficient [ dimregexp, e,-l] - 
al Coefficient[ opregexp, e,-l] ; 
finitediff = 

Coefficient [ dimregexp, e,0] - 
CoefficientC opregexp, e,0]; 

totaldiff = 

pole2diff + poleldiff + finitediff; 

result = 

Integrate [ totaldiff, {x,0,l},{y,0,l}] 

If needed (for speed), the problem can be simplified further by first showing the 
result is independent of p (by showing the derivative with respect to is zero) 
then setting p to zero before proceeding. 
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